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Abstract 

Nonlinear  flow  phenomena  In  transonic  wall  jets 
prototyplc  of  propulsive  lift  devices  such  as 
lifting  ejector  augmenters  and  upper  surface  blown 
wings  have  been  studied  using  the  Rarman-Guderley 
model.  From  modem  line  relaxation  methods,  an 
efficient  computational  method  has  been  developed 
to  treat  the  diversity  of  shock  patterns  produced 
by  various  wall  shapes  and  exit  conditions.  Asso- 
ciated with  the  algorithm  la  a far  field  determlnea 
analytically  from  the  boundary  value  problem  appro— 
P^iete  to  subsonic  conditions  far  downstream. 
Numerical  results  for  circular  arc  boattalls  Indi- 
cate rapid  relaxation  of  the  wall  Induced  disturb- 
ances, even  in  the  supersonic  region.  Partially 
subsonic  and  supersonic  jet  exit  conditions  lead  to 
the  anticipated  wave  Interactions.  Studies  of  other 
shepes  show  that  branch  point  singular  behavior 
associated  with  satisfaction  of  a Kutta  condition  at 
the  wall  trailing  edge  Is  obtained  by  demanding 
continuity  of  the  perturbation  potential  at  this 
point. 

Introduction 

Increased  emphasis  on  propulsive  lift  devices 
in  tactical  and  advanced  aircraft  has  stressed  the 
need  for  greater  understanding  of  the  underlying  ‘ 
fluid  dynamic  processes  which  control  the  degree  of 
force  augmentation  that  can  be  achieved.  Such  con- 
cepts are  exemplified  in  recently  proposed  super- 
critical jet  flap  Implementations  for  advanced 
highly  maneuverable  aircraft  such  as  the  NASA  HlMAT. 
Similar  mechanisms  are  illustrated  in  upper  surface 
blowing  configurations  and  lifting  ejector  aug- 
menters which  are  embodied  In  the  Navy's  XFV-12A.^ 

An  essential  element  in  the  operation  of  these 
devices  are  Coanda-type  wall  jets  consisting  of  jets 
bounded  by  curved  walls  In  which  a transonic  pri- 
mary flow  entrains  an  ambient  secondary  stream 
through  turbulent  mixing  processes.  Existing  models 
for  such  wall  jets  stress  the  Incompressible  treat- 
ment of  these  phenomena  using  eddy  viscosity  and 
energy  methods.  Correspondingly,  there  Is  a need 
for  simulations  that  Include  the  effects  of  non- 
linearities,  mixed  flow,  and  wave  Interactions  on 
the  development  of  the  wall  pressure  distributions 
and  overall  augmentation  forces. 


the  far  field  Is  employed  to  condition  the  numerical 
problem  and  provide  useful  Information  about  the 
decay  of  disturbances.  Both  free  and  wall  jets  are 
discussed  for  several  examples  Illustrating  various 
features  of  this  class  of  flows. 

Formulation 


^^•rtlng  to  the  physical  configuration 
depleted  In  Fig.  1,  a jet  Is  shown  exhausting  from 
the  exit  OC  bounded  by  the  wall  OQ  and  a mixing 
layer  which  has  been  Idealized  as  the  slip  line  CB. 
This  approximation  neglects  turbulent  diffusion 
processes  In  the  study  of  wave  Interactions  with 
the  shear  layer,  but  these  phenomena  can  be 
Incorporated , In  later  refinements.  Furthermore,  It 
will  be  assumed  that  wall  and  jet  turning  angles 
are  small.  In  contrast  to  the  usual  jet  formula- 
tions, in  which  an  upstream  cowl  shape  la  specified, 
or  stagnation  conditions  are  assumed,^"”  this 
analysis  will  treat  a specified  exit  Mach  number 
distribution.  Additional  assumptions  are  Irrota- 
tlonallty  and  subsonic  conditions  Infinitely  far 
downstream.  The  methods  applied  here  can  be 
generalized  to  cases  where  these  restrictions  are 
not  present.  Finite  length  walls  OQ  are  considered 
In  keeping  with  relevance  to  upper  surface  blown 
wings  and  other  propulsive  lift  devices. 

Returning  to  Fig.  1,  the  equations  of  slip 
lines  S.(x,y),  S-Cx.y),  and  the  wall  boundary  B(x,y) 
are  assumed  as 


CB:  Sj^  - y-d  - 6Gj^(x)  - 0 

QA:  $2  " y - iG^(x)  - 0 

OQ:  B - y - 6f(x)  - 0 


where  5 Is  a characteristic  flow  deflection 


a small  disturbance  .limit  in  which 


parameter.  In 

the  scaled  jet  exit  height  D = dd-^'-’,  the  wall  length 
length  L.  ud  the  transonic  similarity  parameter 
K • (1  - are  held  fixed,  as  5 *•  0,  the 

asymptotic  expansions  of  the  velocltv,  pressure  P, 
and  density  p are 


Previous  Investigations  of  related  phenomena 
are  limited  to  the  treatment  of  Invlscid  shockless 
free  jets,  and  Include  the  work  of  Chaplygin. ^ 
Frankl,^  and  Guderley^  all  of  which  employ  hodo- 
graph  methods.  To  study  shock  development  and 
■sixed  flow  phenosMna,  we  have  applied  modem  relaxa- 
tion methods  to  treat  arbitrary  jet  exit  velocity 
distributions  and  assess  the  Influence  of  an  adja- 
cent wall  boundary. 

In  this  paper,  the  computational  model  will  be 
discussed  from  analytical  and  numerical  viewpoints. 
In  analogy  to  unbounded  cases  such  as  airfoil  flows. 


« 


I (x.y;M^i,d.L)  * [l  ♦ «^^\(x,T;K.D,L) 

+ ...jt  + + ...j  j (1.1a) 

P/P»*  1 - + ...  a.lb) 

o/o^*  1 + * ...  (1.1c) 

where  the  subscript  • signifies  conditions  at  x • *, 
♦ is  the  perturbation,potential,  P^  Is  the  ambient 
pressure,  U - a^  - yfjo^  , 
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!•  th«  density,  Sgs  Is  the  speed  of  sound,  q Is  .(he 
flow  velocity  and  the  scaled  coordinate  Y •>  5^'  y 
Is  also  fixed  In  the  limit.  If  a further  transfor- 
mation Is  Introduced  In  which  y ■ Y,  a boundary 
value  problem  can  be  formulated  fot  the  case  of  an 
elliptic  far  field.  The  region  can  be  considered 
as  -the  rectangular  domain  shown  In  ?lg.  2, 
corresponding  to  transfers  of  the  boundary  condi- 
tions to  the  appropriate  undisturbed  streamlines 
allowed  by  the  small  disturbance  limit.  Dropping 
the  tildes,  the  follo%flng  small  disturbance  equa- 
tion holds  Inside  OQABCO 


do  - 6(P,Q)  5 6(x-e)6(y-n) 


where  Q(C,n)  Is  the  source  point  and  P Is  the  field 
point.  Green's  theorem  applied  to  the  region  OQABCO 
gives  the  following  Integrodlfferentlal  equation  for 
♦ 


4 

<(i  - El 


1-1 


1 


where 


5 (d^/dx^  + 3^/3y^)<> 

- (Y^-lJ  (3u^/3x)/2K  , (u  = 4.  ) 


(1.2) 


OO 

Jq  (3u^35)dn 

h"  ♦(C.O)  U (x,y;e.0)dC 


where  we  define  functions  u and  4 such  that 

u(s.y)  " u(x,Y) 

♦(x.y)  - ♦(x,T)  . 


1^  (x,y;e.O)de 


h(n)G(x,y;0,Ti)dn 


Invoking  continuity  of  pressure  and  flow  tangency 
along  the  slip  lines,  we  have,  with  a = 


The  quentlty  G may  be  obtained  from  the  cosine 
transform  where 


G cospxdx 


tlon. 


d(x,a)  - 0 

(1.3a) 

•'o 

♦„(x,a)  - g'.  (x)  , 0 < X 

(1.3b) 

2 /*^  - 

d(x,0)  - Cj^  , I.  < X < «• 

(1.3c) 

G • ^ 1 G cotpxdp 

•'0 

dy(x,0)  - 6*2 (x)  , t < X <“ 

(1.3d) 

Th«  subsldlAry  tqiiatlotM  for 

1 the  constant  C,  la  to  be  computed  by  Itera- 

d^c  2- 

— y - p^G  - 4(y-n)coapC 

In  this  approximation,  tha  slip  lines 

dy 

UWI.  A&Wg  uuluavrwu  AUUWbXWU*  V. 

being  computed  from  the  solution  by  a slsq>le  differ- 
entiation. The  remaining  boundary  conditions  are 


♦y(x,0)  - /(x)  s P(x) 


0 < X < L 


'>,(0,y) 


h(y)  - H(Y) 


(1.3e) 

(1.3f) 


cospC 


t5l. 


'y-n 
5(p,0;5,n) 


y-h’’  y-n- 


C(p,a;^,n)  - 0 


Equation  (1.3f)  Is  representative. of  the  Initial 
a]^t  velocity  profile  which  conceivably  la  deter- 
mined by  the  upstream  duct  contour  and  stvgnatlon 
pressure. 

Far  Field 

To  coeplate  the  formulation  of  the  problem  for 
subaonlc  conditions  far  downstream,  the  asymptotic 
behavior  Is  derived  In  this  section.  Introducing 
a Green's  function  6 satisfying  homogeneous 
Dlrlchlat  conditions  on  OA  and  CB,  and  homogeneous 
tioumana  conditions  on  OC  and  AB  In  Fig.  2 with 


Implying  that 


g . slnhpto}^^Py  cosp^  ^ 

, aiafiBn-aiateltilwsEi  ^ 


Equations  (1.4)  can  be  Inverted  by  a treatment 
of  appropriate  contour  versions  for  the  Inversion 
Integrals,  initially,  without  the  eospC  factor  smd, 
subsequently.  Including  it  using  the  shift  theorem. 
To  obtain  convergence  and  exponential  decay  of  tha 
Integrand,  the  appropriate  cloaure  for  tha  contour 
la  a large  samlelrele  |p|  * K,  I **  ■*,  with  Im  p ^ 0 
for  X ^ 0.  Summing  the  resldMS  at  tha  polea 
p - nvl,  n ■ (agnx)(l,2,3,...),  gives  tha  following 
final  m^reaalon  for  C< 


-irC  ■ •In 

+ ,-no|x-C|| 


a s ir/a 


Equation  (1.5)  la  valid  for  y ^ n.  x % C,  and  can 
bn  nu^ad  aa  cha  Ineagral  of  a (oosacrlc  aarlaa 
giving 


-2TrG  - s(oi(*4<),a(n+y))  + s(a(xK)  .a(n-y)) 

- s(a|x-5|,a(n»y))  - 3(a|x-C|  a(Ti-y)) 

• G - ^ in  ir^ogAT-y)aacha(»4^)-l 

4x  (L^~nn*a(n^)a«cha(rK)  J 


X ri-coaa(n-y)aacha[x-g 
1.1-coaa  (n+y  ) •acha  x-C 


S(A,B)  = £n  ^a”“*coanB  “ -tn|l-a"^| 


j tn(l-2a"^coaB  ♦ a"^) 


Z 5 A + iB 


Aa  Inapaction  of  thaaa  fonulaa  ravaala  that  C is 
axponanclally  saall  aa  x -►  »,  and  is  logarithaically 
singular  at  cha  sourca  point  as  aatlclpacad. 

Basad  on  (1.5),  tha  donlnant  earn  of  tha 
asyapcotic  axpanslon  of  I.  as  x is  glvan  by  tha 
foraula 

* ^^  •"®*sloBy  I slnondn/'  u^alnhci^dC 

Jo  Jo 


vhara  in  tha  avaluaclon,*  cha  contribution  of  Cha 
•iapla  pola  of  vaalshas.  and  intagrala  of  tha 

f Arm  ^ 


A ^ 


f(C.n)dC  and 


.-ci(x+C) 


f«.n)dC 


arlaa.  Tha  ■ultlplicacion  by  u of  cha  aayi^totie 
axpaaalon  raprasancad  by  (1.5)  aa  x -^oand  its 
adbaaquanc  iacagraclon  with  raspaet  to  ( formlly 

* ' 

Tha  uppar  llwlc  of  cha  Innar  iatagral  can  ba  iacar- 
aa  ••  to  vlthla  toraa  of  hlghar  ordar 
lawolvlag  a-«*  aa  x -*  ••.  Thla  Incarpratatloa  la 
nada  in  Iq.  (1.8). 


gives  a development  dominated  by  these  Integrals. 
Writing  the  inner  Integral  of  the  first  of  (1.7), 

r e'®l*"^lf(5,n)d5  - e"®*/*  •®^f(C.n)dC 
•'0  Jq 

+ e“*  f e'®^f(5,n)d«  , 

•'x 

and  if  the  e“®*  factor  of  cha  first  integral  on  the 
left-hand  side  is  Indlcaclva  of  the  behavior  of  p 
••  bban  u*  and  f are  0(a”"*C)  in  this  lladt. 

If  u and  f are  bounded  on  cha  range  of  integration, 
^ |lrac  integral  converges  and  the  second  is 
0(a"20X)  as  x -*  ••.  Tha  second  integral  in  (1.7) 
requires  no  such  dacompostclon  and  is  also  convar- 
gant;  hence,  (l.d)  follows.  Evaluating  cha 
raiwinlng  integrals,  using  (1.5)  and  (l.«,  the 
final  expression  for  the  far-field  is 

* * *FF  " + Cyye"**sinitY* 


0(e”^**)  as  X -*  "• 


(1.8a) 


Sf  ■ ^ f ♦^(5,T)slnhC*dC 

Jo  Jq  * 

- ^ f H(Y)sln7rY*dY  +-^  ♦(C.O). 

•'0  DvT  Jo 


2^  /* 

- ~ I H(Y)sln7rY*dY  +- 
Jo  I 


0)cosh(*dC 


“ 2C^slnh  L* 


(1.8b) 


X*  = irx/Dvf  , C*  5 irC/D^  , L*  5 VL/O^  , 

T*  5 T/D  . 

Humarical  Analysis 

The  nuaarlcal  procedure  is  sisdlar  to  that  one 
first  developed  by  Murman'  and  axcandad  by  JaMSon^ 
and  Bailey  and  Ballhaus.’  Briefly,  the  transonic 
potential  equation  in  divergence  for*  is  diacracisad 
ualng  central  dlffarancas  whan  tha  aquation  la 
elliptic  and  backward  differences  when  it  is  hyper- 
bolic. Thus,  we  nay  wrlta 

(*♦*  - ^ ♦ {*r\ 

• (l-«l^(*ia+l/2-*xl-l/2)  - ^<*Jl+l/2“*rt-l/2^«»l 

‘'l-I^^*xi-l/2"*xl-3/2^  ■ ^^*xl-l/2"*xi-3/2^Pl 
^*Tj+l/2  ■ *T3-1/2^^^J 
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^ L ^ \Pi+l/2  '^1-1/2  /J 

r*  1 

t-*-!  1 i 1-1 
l-’’l+l/2‘*l  ’’l’’i-l/2-* 

u L - ^ WV2\1  • 

^\Pl-l/2  ‘*1-3/2  /J 

rv<-i  Ci-h-z] , < 

L'*l-l/2‘’t  ‘*l‘*l-3/2  J **J 


-3/2 

1-H  *1-1 


'j+1/2  ‘'j’j-1/2 


(2,2) 


*1  ■ ^^‘’l‘’l-l/2  ’ *3  " ^^‘’l‘’i+l/2  • 


and 


•2  - *1  * *3 


Then,  Che  elliptic  difference  expression  In  Eq.  (2.2) 
given  by 


. i 1-1 


,‘’l+l/2^1  ‘‘l’’l-l/2  J 


•1*1-1  ■ *2*1  *3*1+1 


K 


where  Che  missing  subscript  Is  J when  only  I's  are 
present  and  1 when  only  J'a  are  present.  For 
example,  j and  i Also, 


Is  replaced  by 


‘*1  “ ^*1+1  ” *1-1^^^ 


‘‘j  “ ^J+1  ■ ^J-1^^^ 


**1-3/2  " *1-1  ■ *1-2  '‘j-1/2  " ■ ’^j-l 

**1-1/2  “ *1  ■ *1-1 


‘*j+l/2  ■ ^1+1  " 


**1+1/2  " *1+1  ■ *1 


and 


-i: 


If  the  point  (Xj^,Yj)  Is  elliptic 
If  the  point  (Xj^,Yj)  Is  hyperbolic 


Define 


VC^  i K 


. I±i,  /*1+1  ~ *1  ^ *1  ~ '^l-l\ 
^ ' **1+1/2  **1-1/2  / 


Chen 


(0  If  VCj^  > 0 
1 1 if  VCj^  < 0 . 


Her%  the  Iterations  are  viewed  as  steps  In  pseudo- 
time with  4'*'  (NEW)  and  4 (old)  values.  In  addi- 
tion. 


*1  - j+/  1 ^ 1 \ _ ^ / 1 \ 

^1-1/2**!  ^'**1**  1-1/2  **!**  1-3/2'  ^PiP  1-3/2/ 


These  definitions  guarantee  chat  the  linearized 
difference  algorithm  satisfies  the  von  Neumann 
stability  criterion.  See  Jameson^  for  the  proof. 

Ovorrelaxaclon  la  employed  In  the  elliptic 


region  (u 


i.l 


j ■ only.  First,  define 


'l*t-l  - *1*2'“  - *3*1+1 


where  u Is  the  overrelaxaclon  parameter;  l.e., 

1 < 0)  < 2.  Note  that  If  u « 1,  there  Is  no  change 
between  the  two  expressions. 

To  improve  stability  near  the  sonic  points, 
especially  If  a discontinuous  wall  boundary  condi- 
tion Is  being  considered.  It  was  found  necessary,  as 
In  Bailey  and  Ballhaus,^  to  add  to  Eq.  (2.1)  the 
term 


e^t 


xt 


<*t  - *1>  - <*t-l  - *i-l> 

(*1  - x^_^)2 


where  e Is  chosen  to  be  In  Che  range  0 < c ^ .5. 
Boundary  Conditions 

The  boundary  conditions  4(x,D)  ■ 0 and 
4^(x,0)  - f'(x)  may  be  Incorporated  Into  the  numeri- 
cal scheme  using  the  same  techniques  described  In 
Murman  and  Cole.^^  First,  we  shall  concentrate  on 
discussing  Che  boundary  condition  at  Che  jet  exit 
X • 0,  which  may  be  of  one  of  two  types:  (A)  sub- 
sonic at  Che  Jet  exit,  and  (b)  partially  or  com- 
pletely supersonic  at  the  Jet  exit.  The  far  field 
boundary  condition  will  be  created  later. 

(A)  Subsonic  at  the  Jet  Exit 

Here  Che  boundary  condition  Is 


4^(0,Y)  - H(Y) 


0 < Y < D 


where 


K - (Y+1)H(Y)  >0  for  0 < Y < D 


In  this  case,  we  let  Xj^yj  - 0,  Xj^  - 
U-  ■ 0,  and  u,  ■ 0,  and  we  require  x,  - x,  ■ Ax. 
Tnen,  the  x dnlvaclves  In  Eqs.  (2.ljr  and'‘’(2.2) 
becoM 


4 


["^^3/2  - *xl/2>  ■ ^ <«L/2  - ♦L/2>]A» 

• ^ ‘^3/2  ^ ♦xl/2)][*x3/2  - ♦xl/2]/^ 

■ [*  - * «»>)][^  - H(v)]/ax 

(•)  >«rtlxlly  or  Coapl«t«lv  Superxonlc  at  th» 
J»t  Exit 


For  Chia  eas4k  Cvo  boundary  condltlona  are 
required  at  x • 0;  naaaly. 


♦jj(0,Y)  - HCf)  0 < Y < D 


'('^(L-.o)  - ij^a+.o)  - 0 


ia  aatlafled  by  the  aolutlon  for  trailing  edge 
nelghborhooda  In  unmlxed  flow.  Becauae  of  (3),  the 
nonlinear  tern  In  (1.2)  can  be  aaauoed  negligible, 
and  le  locally  hamonlc  In  the  acaled  varlablea. 
Let 

* - X*  + iy  , X*  - X - L 


9 = arg  a , r s mod  a 


w(a)  • u(x*,y)  - lv(x*,y)  ■ coaplax  velocity 


♦ (O.Y)  - g(Y)  for  Y £ [O.D] 


Then  If  the  boundary  condltlona  arc  locally 
llnaariaed  near  the  point  a - 0,  we  obtain 


K - (Y+1)H(Y)  < 0 . 


For  all  polnta  Y £ (O.D)  in  which 
F ” (Y+1)H(Y)  > 0,  g(Y)  need  not  exlat. 

If  a point  (X..Y  ) la  elliptic,  we  uae  the 
dlacratlaatlon  glvCn  in  caae  (A).  On  the  other 
hand.  If  {x-,Y.)  la  hyperbolic,  we  aaatnae  the  grid 
may  be  extefidad  to  the  left  by  Ax/2,  and  we  let 


v(x*,0)  -*  Fa)  = u , X*  t 0 


*0  " *-l  “ •"<*  *03 

Being  Taylor* a theorea. 


♦(0,Y^) 


♦-l,j  ■ *03  - ^♦,(0.Yj)/2  - g(Yj)  - AxH(Yj)/2 

i 

Theae  valuaa  for  ♦.  and  ♦ , , aay  now  be  aubatl- 
tuted  Into  Eq.  (2.23  In  the^dJraal  way,  and  line 
relaxation  aay  be  applied  to  the  flrat  coluan  of 
unknowna  along  x • x^. 

The  far  field  boundary  condition  given  by  Eq. 
(1.9)  contalna  two  unknown  conatanta,  C,  and  C„, 
which  Miat  be  detetalnad  la  an  Itaratlvi  faahlSft. 


u(x*,0)  • 0 


, X*  > 0 


To  doainant  order,  a aufflclent  condition  to  aatla- 
fy  (3)  and  (4)  near  the  origin  la  that  w haa  the 
following  branch  point  behavior 


Hill  + B\/z)  aa  z 0 


(0  < 0 < IT) 


where  B la  a real  conatant  to  be  detemined  by 
matching  with  the  outer  nonlinear  solution.  Equa- 
tion (5)  lapllaa  that 


2 5/2 

♦ - * -uiy  - y Br-"  aln3e/2  (6) 

Several  exaaplas  to  be  dlacuaaed  Indicate  that  the 
approach  deacrlbed  prevloualy  In  which  9 la  maln- 
tainad  contlnuoua  at  the  well  trailing  edge  glvaa 
nuMrlcal  aolutlona  that  aatiafy  the  Kutta  condition 
(3).  However,  a rlgoroua  proof  that  thla  la  an 


The  baale  technique  holda  C_  fixed  while  C,  changea  l»Pllc«tlon  of  the  algorlthe  haa  not  bean  attempted, 
until  the  aolutlon  convergea;'  Than,  C_  la^updated  * alellar  procedure  haa  been  uaad  by  Kruppii  to  aat- 
by  evaluating  the  Integrala  In  Eq.  (inS),  and  the  Kutta  condition  In  the  aolutlon  of  the 

BVdl^AstiaemA  4a  a .e_.  * _ ASIA  1 1 j44  a»<*aWaaaa  v e 


procedure  la  repeated  tmtll  alao  convergea.  In 
order  to  datenalne  C.,  the  eaih  network  la  awept 
froe  left  to  right  ullng  line  relaxation.  Aftar 
the  potential  la  conputad  on  the  line  x - L, 
extrapolation  of  the  Interior  polnta  ylalda 
♦ a,0)  2 ♦ . Aaaunlng  only  that  ♦ la  contlnuoua  at 
X - L,  we  «et  C.  - 9(1,0)  - 9(x,0)  for  x > L,  which 
guaranteea  that  9 (x,0)  - 0 for  x > L. 


It  la  phyalcally  plaualbla  that  a Kutta  condi- 
tion given  by 


lafy  tho  Kutta  condition  in  the  aolutlon  of  the 
tranaonlc  anall  dlaturbanca  lifting  airfoil  problea. 

Reaulta  and  Dlaeuaalon 

In  addition  to  the  aaauaptiona  given  In  the 
Introduction,  the  analyala  prevloualy  deacrlbed  la 
not  directly  applicable  to  choked  flowa  where  up- 
atreaa  and  domatream  condltlona  are  decoupled. 

Sonic  tonea  coeprlalng  the  entire  vertical  dlamnalon 
of  the  flow  field  are  thereby  excluded.  However,  the 
foregoing  methoda  can  be  extended  to  handle  auch 
caaea. 
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A nua^er  of  exaaq)laa  will  now  be  conaldered. 
For  theae  caaea,  the  aaaoclated  wall  dlaplacenant 


functions,  and  exit  velocity  distributions  H(Y)  are 
given  In  Table  1.  For  these  cases,  K is  unity  and 
D will  take  on  this  value  for  the  remainder  of  this 
paper. 


Incompressible  flow  problems  In  which  a singularity 
Is  reflected  between  free  pressure  boundaries 
yielding  an  Image  development  In  which  the  strengths 
alternate  in  sign  to  satisfy  the  slip  line 


Table  1.  Wall  Jet  Cases 


Case 

f ■ Wall  Shape  Function 

H(Y) 

Remarks 

1 

f-fj^-0  ,0<xCl 

- -d-x)^  , 1 < X < 2 

• 3-2x  , 2 < X < L 

.075 

0 < Y < 1 

L > 3.8  for  this  and  all 
other  cases 

2 

f - f^  - -f, 

3 

f-fj-0  ,0<x<l 

- Ix^-l+2(L+l)(l-x)]/2L  , 1<  X < L 

Has  discontinuous  slope 
at  X ■ 1 

■ 

f-f^-0  ,0<x<l 

Has  reflex  curvature  on 
curved  ramp  portion 

5 

f - fj  - -x^  , 0 < X < 1 

- l-2x  , 1 < X < 2 

- -3  , 2 < X < L 

h - 1 0 < Y < 1/2 

- 0 , 1/2  < Y < 1 

♦ • 0.2  on  0 < Y < 1/2 

■ — 

Hall  pressure  distributions  for  the  convex 
rasp  comprising  Case  1 are  shown  In  Fig.  3 for 
K - 1 and  K ■ 1.46.  Interpreting  these  results  as 
those  for  different  final  >^'a  downstream  but  with 
the  same  6,  the  decrease  In  leads  to  upstream 
motion  of  the  terminating  shock  but  a preservation 
of  the  shock  strength.  There  Is  a smooth  accelera- 
tion to  critical  conditions  with  the  location  of 
the  sonic  line  established  at  the  curvature  dis- 
continuity, X ■ 1.  These  calculations,  which  are 
typical  of  the  other  cases,  cost  approximately  $60 
on  the  Berkeley  7600  and  ran  15-30  CP  seconds. 
Figure  4 shows  a close-up  of  the  pressures  near  the 
trailing  adn.  The  dashed  line  has  a slope  propor- 
tional to  *1?  appropriate  to  the  singular  behavior 
given  by  Eqs.  (5)  and  (6)  and  the  Kutta  condition 
(3).  Isobars  shown  In  Fig.  5 are  consistent  with 
these  remarks,  and  demonstrate  the  satisfaction  of 
homogeneous  pressure  boundary  connltions  on  the 
slip  lines.  Because  of  the  weakness  of  the  slngu- 
Isrlty,  e.g. , 0 ~ as  r 0,  special  nuiseri- 

csl  treatments  such  as  those  of  Woods^^  were  not 
used. 

In  Figs.  6 and  7,  rapid  decay  of  the  dis- 
turbances Is  Indicated.  The  relaxation  length  for 
this  decay  from  (1.8b)  with  D " K - 1,  Is  ir.  This 
exponential  decay  la  typical  of  flows  confined  by 
jet  boundaries  and  la  much  more  potent  than  for 
bodies  In  unbounded  fields.  Qualitatively  similar 
effects  have  bera  discussed  by  Murman,^^  and 
Plntola  and  Lo^  In  connection  with  tunnel  wall 
Inter ferenca  un  transonic  airfoils.  The  distinc- 
tion between  confined  and  unconflnad  flows  can  be 
appreciated  by  an  Intarpretaelon  of  the  exponential 
eerlas  dua  to  (1.3)  arlalng  In  the  far  field 
davalopmants  dominated  by  (1.8).  Because  of 
linearity  of  the  far  flel4  this  sarlea  la  directly 
related  to  expansions  occurring  In  analogous 


condition.  Thus,  the  relaxation  to  uniform  condi- 
tions downstream  which  must  be  consistent  with  homo- 
geneous conditions  on  the  slip  lines  produces  a more 
rapid  decay  than  found  In  unconflned  flows. 

In  Figs.  8a  and  8b,  the  upper  and  lower  slip 
lines  obtained  from  Integration  of  (1.3b)  and  (1.3d) 
are  given  for  Case  1.  It  is  evident  from  Fig.  8a 
chat  Che  curved  surface  In  this  approximation  turns 
the  flow  so  that  the  streams  are  parallel  for  x "*  *•. 
In  the  near  field.  It  Is  evident  from  Fig.  8b  that 
this  Is  not  quite  the  case.  Asymptotic  parallelism 
can  be  established  for  subsonic  conditions  far 
downstream  by  Integration  of  the  small  disturbance 
continuity  equation.  Thus 

g'^(x)  - ^ (Ku(x,y)  - ^ u^)dY  (7) 

and  since  u ■ 0(e”**) , this  expi'esslon  shows  that 
GjCx)  ■•'CjCx)  as  X Equation  (7)  can  also  be 

obtained  by  differentiating  (1.8a)  with  respect  to 
Y and  using  (1.3b)  and  (1.3d).  For  the  case  of  a 
free  jet  with  a symmetrical  exit  Mach  number  profile 
function  H,  l.e.,  H(Y-l/2)  - H(l/2-Y),  Cj^  - -C^  and 
the  divergence  theorem  or  Integration  of  (7)  between 
X ■ 0 and  " gives  the  displacement  of  the  Jet 
Infinitely  downstream  as 

Cl«  - 7 [kH(Y)  H^jdT  (8) 


• 

Such  a series  can  be  suaeed  by  recognising  that  it 
la  a partial  fraction  expansion  of  a hyperbolic 
function  which  Is  exponentially  small  as  x **  **, 
consistent  with  developments  such  as  (1.8a). 
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analytical  solution  obtained  either  from  sunsned 
eigenfunction  expansions  or  transfonsc  la: 


where  0^(0)  - 0 has  been  used. 

Insight  Into  the  SMchanlsiss  causing  accelera- 
tion to  supercrltlcallty  can  be  obtained  from  the 
flow  direction  field  for  Case  1 shown  In  Fig.  9. 

For  clarity » all  Isocline  slopes  have  been  magni- 
fied a factor  of  100,  and  only  the  entrance  sec- 
tion 0 < X < 1 Is  depicted.  The  expansion  around 
the  curved  ramp  on  the  interval  1 < x < 2,  leads  to 
upstream  Influence  In  the  subsonic  region  which 
turns  most  of  the  flow  downward  upstream  of  x ■ 1, 
producing  throats  and  acceleration  due  to  the 
stream  tube  contraction  required  by  the  zero  slope 
boundary  condition  In  that  region.  Also  Indicated 
f*  the  "ballooning”  due  to  the  singularity 
occurring  at  (0,1),  the  top  point  of  the  Jet  exit 
station.  In  contradistinction  to  the  trailing  edge 
where  u -•  0,  local  linearization  cannot  be  used  to 
characterize  the  flow  behavior  In  this  region,  and 
some  local  similarity  solutlra  must  be  sought, 
presumably  of  the  form  d - x“f(y/x8)  where  a and  B 
are  exponents  to  be  determined. 

By  contrast,  the  concave  shape  shown  In  Fig. 

10  produces  the  anticipated  compressive  decelera- 
tion, which  Is  also  depicted  In  Fig,  11, 

The  effect  of  a slope  discontinuity  is  indi- 
csted  In  Fig.  12.  It  Is  evident  that  the  nusari- 
cal  method  accurately  locates  the  initiation  of  the 
sonic  line  at  the  point  (1,0)  where  the  flow  Is 
’tripped"  to  criticality  by  the  acceleration 
singularity  at  this  location.  In  moat  othsr 
respects,  the  pressure  distribution  is  similar  to 
chat  for  Casa  1. 

I®  Wg.  13,  the  effect  of  a reflex  curvature 
In  accelerating  cha  recompresslon  process  Is  shown. 
As  related  to  comparable  Cuming  and  wall  deflec- 
tion treated  In  Case  1,  the  strength  of  the  ter- 
minating shock  is  considerably  Increased  as  well  as 
the  magnitude  of  cha  pressures  near  the  trailina 
edge.  * 


effect  of  mixed  flow  conditions  at  the 
aheam  In  Figs.  14-16.  Here,  cha  function 
H as  well  as  ®(0,Y),  comprise  the  Cauchy  data 
needed  to  properly  pose  the  hyperbolic  portion  of 
the  Initial  manifold.  Since  the  vertical  velocity 
«y(0,Y)  can  be  derived  as  a tangential  differen- 
tiation, the  Cauchy  data  connotes  specification  of 
Che  additional  velocity  component  for  supersonic 
portions  of  the  Jet  exit  station.  Figure  16  Indl- 
cates  chat  In  addition  to  the  usual  terminating 
shock,  the  transition  from  hyperbolic  to  elliptic 
flow  occurs  across  a weak  shock  emanating  from  tha 
specified  u discontinuity  at  (0,1/2), 

In  Fig.  17,  tha  behavior  of  centerline 
pressures  for  various  free  Jet  cases  Is  shown.  The 
monotone  subcrltlcal  behavior  exhibited  by  these 
nonlinear  cases  has  not  yet  bean  corroborated  by 

proof  based  on  Cha  boundary  value  problem 
with  subsonic  exit  and  downstream  conditions.  For 
linearized  subsonic  flow,  this  property  Is  obvious 
from  Che  maximum  modulus  theorem.  It  should  be 
nbtad  that  In  cha  free  Jet  problem  discussed  here, 
specified  mssa  flow,  pressurs  ratio  and  final  Mach 
number  uniquely  datamlna  d — the  scale  paramscer 
for  Che  Jet  displacement. 

• validation,  a coiq>arlson  of  numerical  and 
^andtl-Clsuert  free  Jet  solutions  for  u(x,l/2)  for 
H - 0.35  and  K - 10  la  shewn  In  Fig.  18,  where  the 


u 


V 


and 


n 


slny^ 

slnhx^ 

coshx^  cosv^ 
coahx  - cosy 


x^  = itx/vic 


(8a) 


(8b) 


The  slight  discrepancies  shown  in  Fig.  18  presumably 
derive  from  the  small  nonlinear  effect  associated 
with  tha  finite  R value,  and  truncation  errors  of 
f^e  discretizations  which  are  only  approximately 
second  order  for  a non-unlfom  elliptic  mesh.  The 
associated  universal  slip  line  curve  Is  obtained 
from  the  following  Integral  of  (8b) 


Cj(x) 

HR 


1 _ [1  ^^■«(2n^»^l)le~^^”^*^^^* 

^ I 0 (2n+l)^ 

- j ctnh”^(coshx)| 


(9) 


where  the  daggers  have  been  dropped.  Equation  (9) 
has  Che  following  limiting  behavior 


Cj(x) 

HR 


xJnx 


ln4-*-l|x  0(x^)| 


as  X -*  0 


Cl  («)  , . _ 

HR  “ T ■ ~ as  X 

IT 

and  is  plotted  in  Fig.  19.  The  asymptotic  half 
width  thus  checks  that  given  by  (8)  whan  the  non- 
linear second  term  In  the  integral  la  omitted. 

ConcluelonH 

Submerged  Invlscld  transonic  wall  jets  have 
been  created  for  subsonic  conditions  Infinitely  far 
downstream  for  a specified  exit  velocity  profile. 
Unchoked  flow  has  been  assumsd  and  a series  of  wall 
shapes  have  been  considered  giving  rise  to  the 
following  observations: 

• Line  relaxation  provides  an  efficient  means  to 
treat  the  diversity  of  wave  patterns  chat  can 
occur,  particularly  for  nixed  Initial  condi- 
tions. Accordingly,  Inverse  probleae  In  which 
contours  are  Identified  Co  reduce  wave  drag  ar1 
enhance  thrust  recovery  nay  represent  a feasi- 
ble near-term  possibility. 

• A Ructa  condition  can  be  satisfied  merely  by 
requiring  continuity  of  the  potential  across 
the  trailing  edge.  The  numerical  solution 
tracks  a local  singular  solution  which  has  a 
square  root  zero  and  Is  locally  harmonic  In 
scaled  variables. 

• As  compared  to  unconflned  flows,  the  slip  line 
boundaries  create  a rapid  decay  of  the  dis- 
turbances. The  functional  form  of  the  far 


field  perturbation  potential  Is  an  exponen- 
tially damped  sine  similar  to  Incompressible 
flow  but  different  In  that  Its  amplitude 
Interacts  nonllnearly  with  the  near  field. 

• Analysis  of  the  free  Jet  case  Indicates  sub- 
crltlcal  monotone  streamulse  variations  of  the 
pressure  for  a subsonic  Jet  exit,  as  In  linear 
subsonic  flow. 

• Acceleration  of  the  wall  jet  to  criticality 
over  convex  walls  Is  accomplished  by  stream 
tube  contractions  and  throats  Induced  by  up- 
stream Influence  of  the  turning. 

• The  upper  slip  line  of  the  wall  Jet  becomes 
asymptotically  parallel  to  that  emanating  from 
the  trailing  edge  Infinitely  far  downstream. 
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Fig.  1 Wall  Jet  configuration. 


Fig.  2 Transformed  wall  Jet  domain. 


Fig.  3 Pressures  along  wall  for  K • 1 and  K • 1.46, 
fj^  ■ normalized  wall  ordinate. 
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6 Cp  surface  for  Case  1 


Fig.  7 Mach  surface  tor  Case  1 


Fig.  8b  Closeup  of  slip  lines  near  Jet  exit  for 

Case  I. 


Pig.  8a  Slip  line  shapes  for  Case  1 


Fig.  15  Pressure  surfsco  for  Case  5 


Fig.  14  Pressures  along  wall  for  Case  5 (partially 
supersonic  jet  exit),  f.  ~ normalized  wall 
ordinate. 


Pig.  16  Mach  contours  for  Case  5 


Centerline  pressures  for  a transonic  free 
jet  for  various  pressure  ratios  and 
similarity  parameters,  K. 
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lead  to  the  anticipated  wave  interactions.  Studies  of  other  shapes  show  that 
branch  point  singular  behavior  associated  with  satisfaction  of  a Kutta  condition 
at  the  wall  trailing  edge  is  obtained  by  demanding  continuity  of  the  perturbation 
potential  at  this  point. 
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